Properties of parallel upper critical field within Continuous Ginzburg-Landau model 
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In this paper, we employ a continuous Ginzburg-Landau model to study the behaviors of the 
parallel upper critical field of an intrinsically-layered superconductor. Near Tc where the order 
parameter is nearly homogeneous, the parallel upper critical field is found to vary as (1 — T/Tc)^^^. 
With a well-localized order parameter, the same field temperature dependence holds over the whole 
temperature range. The profile of the order parameter at the parallel upper critical field may be of 
a Gaussian type, which is consistent with the usual linear Ginzburg-Landau theory. In addition, the 
influences of the unit cell dimension and the average efi'ective masses on the parallel upper critical 
fleld and the associated order parameter are also addressed. 



I. INTRODUCTION 

Most high Tc superconductors (HTSs) have layered 
structures and the layered superconductivity is closely 
related to the behavior of the order parameter. Spatial 
variation of the order parameter would aid us to intu- 
itively understand various properties such as the coher- 
ence lengths of layered superconductors. On the other 
hand, the investigation of the upper critical field Bc2 
may provide information on the coherence lengths, and 
in principle, allow the testing of existing theories (for 
example, see Refs. 0-0). Hence, a study involving the 
order parameter and the upper critical field should prove 
valuable. 

The phenoinenological continuous Ginzburg-Landau 
(CGL) modela is convenient for a description of the or- 
der parameter and the upper critical field of layered su- 
perconductors. The coefficients in the CGL free energy 
are assumed to be spatially dependent. As a result, the 
amplitude of the order parameter varies, reflecting the 
layered nature. The amplitude of the order parame- 
ter at weakly superconducting layers may be extremely 
small, corresponding to a weakly linked layer system 
similar to the Lawrence-Doniach (LD) model. On the 
other hand, when the spatial dependence is neglected, 
the CGL model is reduced to the anisotropic Ginzburg- 
Landau (GL) model. Hence, the CGL model approaches 
the limiting cases of the LD model and the anisotropic 
GL model. ^ 

In a previous workEI, a set of spatial coefficients for 
the CGL model was proposed for a layered supercon- 
ducting system in which the unit cell was assumed to 
compose of equivalently thick superconducting and insu- 
lating lancers and no applied magnetic fleld was present. 
RecentlyO, a magnetic field parallel to the layers with 



different thickness was introduced and efficient comput- 
ing methods have been adopted to determine the generic 
properties of the parallel upper critical field B^^ of var- 
ious layered superconductors. In the present work, we 
shall examine various features pertinent to 3^.2 and the 
associated order parameter of a typically-layered super- 
conductor. 



II. MODEL 

In the CGL model of Koyama et aZ.i, layered supec.- 
conductors have been classified into three categorieaij, 
one of which the layered HTS BiaSraCaCusOg (Bi2212) 
may fall into. Following Ref. ||, Bi2212 shall be chosen 
as our modeling prototype for laired superconductors 
as it possesses a large anisotropyM, and thus is suitable 
for a detailed study that examines the relationship be- 
tween the amplitude of the order parameter and the lay- 
ered structure. Note, however, that a study involving 
the phase effect jof the order parameter ( Josephson cou- 
pling) in Bi2212E3 may require a LD descriptioru for an 
appropriate investigation. 

The unit cell of Bi2212 composes of two Cu02 bilayers, 
separated by the BiO-SrO interlayer, which is referred to 
as insulating (I) layer for convenience. The two adjacent 
Cu02 planes of the bilayer (interplane distance ~ 3 A) 
are strongly coupled so that they can be treated as a 
single superconducting (S) layer; therefore the distance 
between two supercorulucting layers is half the c-axis lat- 
tice constant ~ 15 M3. Denoting the thickness of the I 
and S layers as dj and c?5, respectively, we may write the 
size of the unit cell D ^ D = dj + ds- The CGL free 
energy for the system is 
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with planar vector r = {x,y) and vector potential 
A{f,z) = (A(2)(f, z),A2(r, z)). M(z) denotes the effec- 
tive masses along the z-direction (c-axis), and m(z) is 
the corresponding planar parameter. The GL coefhcient 
a{T, z) and the effective, masses are taken as before (/3 is 
assumed as a constant )l3: 

a{T, z) = [ao + ai cos{2iizlD)] (1 - T/TJ, (2a) 
Go + GiCos(27rz/i:>), (2b) 
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where ao, ai, Go, Gi, go and 171 are the model parame- 
ters. These parameters are found to be related to exper- 
imental data, to each other, and to the intrinsic param- 
eters, di and ds (see discussions in section IV). 

Let an external magnetic field B be applied parallel to 
the a- or &-axis which is in the y-direction. Thus, the 
vector potential can be taken as A — {Bz, 0,0). Assum- 
ing ^'(f, z) = e'''ii ''vl/(z), it follows from Eq. |^ that, by 
minimizing the free energy, 
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with zo — hkx/{2eB). At B = Bc2, the superconducting 
order develops in the S layer first so that one may choose 
Zo = D/2. The high order term in Eq. | I3\^ {z)\'^^ [z) , 
may be omitted since the order parameter at i?c2 is phys- 
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For a given temperature T, the maximum magnetic field 
B which satisfies the above equation, together with the 
boundary conditions *(0) = *(-D) and ^*(-2)|^^£, = 0, 
gives a point on the Bc2-T curve. The largest B can be 
readily ajchieved by treating B^ in Eq. ^ as eigenvalue 
problemaj. 



III. RESULTS AND DISCUSSION 

The order parameter distribution in a unit cell at dif- 
ferent temperatures is plotted in Fig. 1(a). At low tem- 
peratures, the order parameter is mainly confined within 
the S layer, signifying a two-dimensional (2-D) state. At 
high temperatures, it effectively penetrates into the I lay- 
ers. Near Tc, it varies smoothly and is nearly a constant 
throughout the unit cell, indicating a three-dimensional 
(3-D) state behavior. The present model thus correctly 
accounts for the behavior of a 2-D state at lower tem- 
peratures and a 3-D state near Tc. Note that the weak 
modulation of the order parameter near may generate 
a genuine 3-D superconductor. This i&,diffcrcnt from the 
so-called 3-D region of the LD modelEJ, where the inter- 
layer coherence length is much larger than the interlayer 



ically small. To explore the features of the order param- 
eter along the z direction, we assume ky = 0. Finally, we 
obtain 
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spacing or the size of the unit cell, but the order param- 
eter is still assumed to be discontinuous. Thus, there 
is the possibility of a true 3-D superconductivity with a 
nearly uniformly distributed order parameter even in a 
highly anisotropic superconductor. This situation can be 
obtained by just varying the temperature (see Fig. 1(a)). 

Again, it is found that the peaks of the order parame- 
ter can be fitted by a Gaussian function. The exponen- 
tial factor is the most significant part of the Gaussian fit, 
showing that the ground state of the CGL linearjenuation 
is similar to that of the usual linear GL equationoEj. We 
emphasize that this similarity, together with many rea- 
sonable results to be presented, reveals the plausibility 
of our methods of calculating Bc2- The fitted ^±(0) is 
0.96 A, which compares favorably with some experimen- 
tal values of ~ 1 M3. 

The calculated parallel Bc2 as a function of temper- 
ature is shown in Fig. 1(b) and Fig. 2. Near T^, the 
feature of Bc2-T is square-root like while far away from 
Tc, it is linear. The linear behavior in Fig. 1(b) can 
be understood by identifying B^ oc 1/^^ while the lat- 
ter is proportional to 1 — T /T^. Note that the rela- 
tionship between Bc2 and T is also linear within the 
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anisotropy GL theory, in which B'^^i^) = 27r^ {t)^±(t) ' 
where both the interlayer and in-plane coherence lengths 
^±{T) and C||(T) are proportional to (1 - T/Ta)-^^^ so 

that bI^{T) oc (1 -T/Tc). 

Since the parallel i?c2 in Bi2212 rapidly exceeds acces- 
sible laboratory magnetic fields when the temperature is 
reduced from Tc, only the c alcula ted data near Tc can be 
compared with experimentstZltS (see Fig. 2). By consid- 
ering a constant solution of the order parameter to Eq. 4, 
one can immediately obtain a square-root Bc2-T relation 
near Tc- Note, however, that with open boundary condi- 
tions (OBC, 'I'(z)|z^±oo = 0) imposed on Eq. 4, we have 
found that there is a linear Bc2-T relation near Tc (see 
solid circles in Fig. 2). The deviation from the linear be- 
havior can be understood as a dimensional crossover (see 
Ref. |l^ and references therein). Moreover, when using 
spatial-independent coefficients (AGL) in the OBC sim- 
ulation, we obtain a linear Bc2-T relation in the whole 
temperature range (as expected). Note that these calcu- 
lation results are in reasonably and qualitatively agree- 
ment with experiments (which are diverse due to factors 
such as crystal quality, measurement methods, etc). 

The square-root behavior near Tc indicates that the 
upward curvature of the Bc2-T curve iS|-absent in the 
present simulated system. Recent studiesHlHS show that 
the upward curvature is perhaps not intrinsic. Indeed, 
such curvaturft-ia neither found nor obvious in the WHH 
approximatiorO, t^ d-wave theorycj and the mixed d- 
and s-wave theoryu. Note, howeverp. that the feature of 
the Brs^T curvature is controversialQ and remains to be 
testednJ. According to Fig. 2, the curvature of Bc2-T is 
boundary-condition dependent and thus indeed difficult 
to arrive at an absolutely conclusive conclusion. Tba 
curvature may also be affected by physical phenomenal 
such as the spin orbit scattering (for example, see the 
microscopic theory in Ref. ||). 

The calculated Bc2 at zero temperature is about 700 
Tesla and is comparable with those extrapolated from eju 
periments on Bi2212B|dand other HTS such as YBCOEI 
and Bi2Sr2Ca2Cu30ioE3. However, the expierimentally 
extrapolated data might not be too reliablenS and a pos- 
sible new way ta detect Bc2 with Josephson plasma has 
been suggestedc3. 

Fig. 3(a) shows the spatial distribution of the order 
parameter for different D of the unit cell at zero tem- 
perature. For small D, the order parameter covers the 
entire I layer but as D increases, its penetrations into 
the neighboring I layers become restricted and quickly 
fall to zero. The 3-D (2-D) behavior for small (large) 
unit cell is in accordance with the features reflected 
in the calculated Bc2-T curves of Refs. ||,^. Fig. 3(b) 
presents the variation of the upper critical field with the 
size of the unit cell at T = K. The obtained criti- 
cal field decreases with_/3, which is qualitatively consis- 
tent with the g3 modeH. Here a power law (dash line) 
could not fit the trend well but an exponential fit (solid 
line) is acceptable. Experimentally, for the similar com- 



pounds of Tl2Ba2Ca„_iCu„02n-i-4 (n=l,2,4), Mukaida 
et alEB reported that the upper critical field generally 
decreases as the number of Cu02 layer increases. They 
attributed their results to the effects caused by the dif- 
ferent thickness of the effective superconducting layers in 
Tl2Ba2Cu06, Tl2Ba2CaCu208 and Tl2Ba2Ca3Cu40i2, 
whose respective c-axis lattice constants are 23.2, 29.3 
and 41.9 A. Clearly, the theoretical trend presented in 
Fig. 3(b) is consistent with their experimental observa- 
tions. 

The mass dependences of the parallel upper critical 
field at zero temperature are shown in Fig. 4(a) and (b). 
In these calculations, the value of Gi (gi) was fixed while 
that of Go (50) was varied to obtain the varying AI± 
(to||). It is clear that large values of both M± and m|| 
result in a large critical field, which is. consistent with the 
anisotropic GL theory, the LD modelEl and the ga theoryO. 
It is worth mentioning that as Mj_ increases, we find that 
the order parameter ^'(D/2) at the S layer grows while 
that at the I layer, \E'(0), decreases, leading to a larger 
difference of - ^'(O). Since l/[^{D/2) - 4'(0)] is 

approximately proportional to the strength of interlayer 
coupling between adjacent S and I layersu, thus M± sup- 
presses interlayer coupling. In contrast, TO|| is found to 
enhance interlayer coupling. We attribute this to the ef- 
fect that M± enhances the order parameter in the super- 
conducting layer while TO|| suppresses it. Hence, HTS are 
intrinsically favorable for a large M±, which corresponds 
to a weakly linked layered system. 

The spatial variation of the order parameter at sev- 
eral condensation energies is plotted in Fig. 5(a). The 
temperature is set to zero. It is obvious that the largest 
energy corresponds to the largest order parameter in the 
S layer and that the smaller the condensed energy, the 
broader the order parameter. Fig. 5(b) further shows 
that the critical field increases with the condensation en- 
ergy 

Up till now, we have set the parameters ag, ai, Gq, 
Gi, go and gi (see Eq. 2) in the S layer the same as those 
in the I layer (for example, ap = cxq)- We shall now con- 
sider the case where Q!q ^ ctg. The ratios of q^q /o^o 
Fig. 6 (a) and (b) are 20 and 2000000, respectively. The 
temperature is zero. The non-monotonic trend (first as- 
cending and then descending) in Fig. 6(a) is qualitatively 
consistent with the data extracted from the numerical 
work of Refs. ||. 

When the energy condensed in the S layer is extremely 
large, the system would be in an extreme 2-D state, as 
illustrated in Fig. 6(b) for different unit cell sizes. The 
order parameter totally resides in the S layer. The S layer 
fully decouples with the adjacent I layers and therefore 
the system is in an extreme 2-D state. This interesting 
2-D behavior can be confirmed by the thickness depen- 
dence of the upper critical field. The theoretical data can 
be fitted by an inverse relation typical of a 2-D system 
(for example, see Refs. ||,|l^). It is interesting to find that 
the 'i'{D/2)-D profile is qualitatively consistent with the 
Bc2-D trend. Such qualitative consistency can also be 
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found in Figs. 6(a) and 3. 

The extreme 2-D behavior can be further confirmed 
by the spatial distribution of the order parameter and 
the temperature dependence of the upper critical field, 
which are shown in Fig. 7. In Fig. 7 (a), the order pa- 
rameter drops down sharply and is confined in the S layer 
in a large temperature range till 84.9 K. In Fig. 7 (b), a 
square-root relation between B1.2 and T holds over the 
whole temperature range and this again is a typical 2-D 
behavior, which has been reported in the literature (for 
example, see Refs. f|,p|J3|j3|). 



IV. CONCLUSION 



Within a continuous Ginzburg-Landau model for lay- 
ered superconductors, we have calculated the parallel up- 
per critical field and the associated order parameter with 
respect to the variation of the temperature, the unit cell 
dimension, the average effective masses and the GL con- 
densation energy. Near the vicinity of Tc where the order 
parameter is nearly homogeneous, the parallel upper crit- 
ical field is found to be square-root like. With a highly 
localized superconductivity, the same field temperature 
dependence holds over the whole temperature range. The 
order parameter at Bc2 of the linear CGL equation may 
demonstrate a Gaussian profile, which is consistent with 
that of the usual linear GL equation. The profile of the 
maximum order parameter in the superconducting layer 
against the unit cell size may be correlated with the trend 
of the upper critical filed versus the unit cell dimension. 
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FIG. 1. (a) Spatial distribution of the order parameter for 
different temperatures, (b) Temperature dependence of the 
parallel upper critical field. 
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FIG. 2. Calculated temperature dependences of the par- 
allel upper critical field near Tc, compared with experiments. 

The solid squares correspond to the periodic boundary condi- 
tion with spatial-dependent coefficients (CGL), the solid cir- 
cles to the open boundary conditions with spatial-dependent 
coefficients (CGL) and the open diamonds to the open bound- 
ary conditions with spatial-independent coefficients (AGL). 
The solid hno is a fit varying as (1 - T/Tc)°-''. The dotted 
line signifies the crossover temperature from 3D to 2D in the 
PBC-CGL (solid squares) and OBC-CGL (soUd circles) cal- 
culations. 

FIG. 3. (a) Order parameter distribution and (b) upper 
critical field at different sizes of the unit cell. 



FIG. 4. Upper critical field for (a) perpendicular average 
mass M± and (b) parallel average mass my. 

FIG. 5. a) Order parameter distribution and (b) upper 
critical field at different condensation energies. 

FIG. 6. Order parameter and upper critical field for 
ai/a'o = 20 in (a) and ai/a'o = 2000000 in (b). The profiles 
of the maximum order parameter vs D seem consistent with 
the corresponding Bc2-D trends. The dash line in (b) is ap- 
proximately an inverse fit while the solid line is an exponential 
decay. 

FIG. 7. Typical 2D temperature dependence of the paral- 
lel upper critical field. 
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